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Abstract 

We show that fluctuating proton positions in the colliding nuclei generate, 
on the event-by-event basis, very strong magnetic and electric fields in the 
direction both parallel and perpendicular to the reaction plane. The mag- 
nitude of E and B fields in each event is of the order of ^ 10^^ Gauss. 
Implications on the observation of electric dipole in heavy ion collisions is 
discussed, and the possibility of measuring the electric conductivity of the 
hot medium is pointed out. 
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1. Introduction 

The importance of the strong magnetic field created in heavy ion collisions 
has been recently emphasized in Refs. [l|, S [sj , where the chiral magnetic ef- 
fect has been proposed and studied. The chiral magnetic effect results in the 
electric current along the direction of the magnetic field B. Consequently, 
the measurable charge separation may be observed in the direction perpen- 
dicular to the reaction plane {y axis in Fig. [p - the dominant direction of 
the magnetic field averaged over events [jj, Recently, the STAR collab- 
oration has published |5| first results on the charge-dependent two-particle 
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correlations with respect to the reaction plane, which are qualitatively con- 
sistent with the chiral magnetic effect expectationsjthis conclusion, however, 
is under intensive debates (see, e.g., Refs. BHiiliol)- 

As was mentioned earlier, according to early estimates (see Refs. [H, 0]) 
the magnetic field perpendicular to the reaction plane in heavy ion collisions 
is dominant and its magnitude may reach very high values up to dozens of 
^ 10^^ Gauss. However, in all previous analysis the fluctuation of the 
magnetic field in each event was disregarded, and instead the averaged value 
over many events was used to calculate (B). The importance of the event-by- 
event calculation becomes apparent after calculating the x component of the 
magnetic field at x = (denoted by a black dot in Fig. [T]). Indeed, averaging 
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Figure 1 : The transverse plane of a peripheral heavy ion collision with the impact param- 
eter b. We calculate the magnetic B and electric E fields at t = (time of the collision) 
and X = (denoted by a black dote) and we also average over a transverse area denoted 
by a grey circle. 

over many events we obtain (Bj.) = (as we explain later), suggesting that 
the X component of the magnetic field can be neglected. In this Letter, we 
argue that owing to fluctuating proton positions in both nuclei, B^ becomes 
comparable to By on the event-by-event basis. We also calculate the electric 
field E and show that the x and y components of B and E fields are of the 
same order of magnitude, at least at the early stage of the collision, when 
the hot medium response may be neglected. 

In the next Section we describe our approach in detail and present our 
results. Comments and summary are presented in Section 3 and Section 4, 
respectively. 
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2. Calculation 



The electric and magnetic fields at a position x and observation time t 



can be calculated in each event according to the following equations [H 



V- 1 - 



eE{t,x) = aEM^ r^^n, 



1 — v"^ 

eB{t,x) = OeM^ -jj^VnXRn, (1) 

" Rl (l - X Vn]yRl) 

where sums are performed over all protons in both nuclei. The fine-structure 
constant oem = e'^/in ~ 1/137 and Rn = x — Xn(t), where Xn is a position of 
proton moving with the velocity Vn- If Vn is non-zero only in the z direction, 
then [Rn x tVi)^ = Rn±v^z^ where the transverse vector between a proton 
and the observation point Rn,i. = x± — Xn,± does not depend on time 

In our calculation we sample proton positions a?^ at t = according to 



the Woods-Saxon distribution with the standard parameters [1^. We also 
assume that both nuclei are infinitely thin; and that all target and projectile 
protons move with the same velocity v^^^^ = [0, 0, Vz] and v^^"^ = [0, 0, —Vz], 
respectively. The value of Vz is defined by the collision energy a/s and the 
proton mass nip, v1 = 1 — {2mp/^/s) . Finally, when calculating B and E 
at a given point x we take into account only those protons which are not 
closer than rcut = 0.3 fm to the observation point x. This cut-off is to be 
implemented owing to the singularities of Eqs. ([1]) at i?„ — )■ 0. The value 
Tcut = 0.3 fm was fixed as an effective distance between partons in a nucleon. 
By a variation of r^ut in the range from 0.3 fm to 0.6 fm we affirm weak 
cut-off dependence of the fields. 

We performed our calculations at t = and at a; = (denoted by a black 
dot in Fig. [1]). The results for magnetic field in AuAu collisions at ^/s = 200 
GeV are presented in Fig. O 

As seen from Eq. (IT]) {B^) = at x = ol However, in each event the x 
component of the magnetic field can be huge and vanishes only after averag- 



-'^It should be noted that the retardation effects are already taken into account and i?„ 
depends explicitly on the observation time t. 

^Indeed, one nucleus gives {B^) ~ (yn) ~ 0, in contrast to (By) ^ (a;„) 7^ if & > 0, 
see Fig. [TJ 
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Figure 2: The mean absolute value of the magnetic field at i = and a; = as a function of 
impact parameter b for AuAu collision at \/s ~ 200 GeV. Fluctuation of proton positions 
lead to non-zero values of (|-Bx|) that are comparable to {\By\). 

ing over events. To study the magnitude of the magnetic field disregarding its 
direction from one event to another we consider the average (over events) ab- 
solute valu^of the magnetic and electric fields, {\Bx,y\) and {\Ex,y\). Due to 
fiuctuations of the proton positions we obtain comparable numbers for (l-Bxl) 
and (I -By I) suggesting that on the event-by-event basis we should expect huge 
fields both in x and y directions. Since the chiral magnetic effect leads to the 
electric current along the magnetic field, our result indicate that, in principle, 
the chiral magnetic effect may take place not only in the y direction but also 
in the x direction. 

The results for the electric field are shown in Fig. [31 

The symmetry of the system presented in Fig. [1] implies that at x = 
the average value of the electric field (E^) = {By) = 0. However, as seen in 
Fig. [31 fiuctuations lead to (l-Exl) ~ {\Ey\) with magnitude of the order of 
m^. It is interesting to notice that x and y components of the electric field 



■^First we calculate, e.g., in an event from all protons and after that we take the 
absolute value. Next we calculate average over events. 
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Figure 3: The mean absolute value of the electric field at t = and a; = as a function of 
impact parameter b for AuAu collision at y/s = 200 GeV. Fluctuation of proton positions 
lead to non-zero values of (|i?a:|) ~ {\Ey\). 



are almost identical to the x component of the magnetic field 

{\B^\) ^ {\E^\) ^ {\Ey\) . (2) 

This conclusion can be drawn easily from the analysis of Eqs. ([1]). The fields 
Bx and E^^y are driven only by fiuctuations of the proton positions and since 
Vz ^ 1 it is evident that Eq. (jj]) is valid. For peripheral collisions, the y 
component of the magnetic field is infiuenced not only by the fiuctuations, 
but also (mainly) by the geometry of the collision, as seen from Fig. [21 where 
{\By\) and (By) are compared. Thus we expect {\By\) to be larger than the 
fields in Eq. 

To provide more complete information on the fiuctuating values of B and 
E fields, we show in Fig. IHthe event-by-event histograms for B^^y and E^^y 
at the impact parameter 6 = 8 fm. As expected By distribution is shifted 
away from zero and B^, E^ and Ey distributions are practically indistin- 
guishable, consistent with Eq. ([2]). We checked that for 5 = histograms 
look very similar, the only difference is By that is centered around zero being 
indistinguishable from B^, E^ and Ey. 
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6 = 8 /m, AuAu, VJ = 200 GeV 




Figure 4: Normalized event-by-event histograms of x and y components of magnetic B 
and electric E fields at t = and af = at the impact parameter 6 = 8 fm. As expected for 
peripheral collisions, By distribution is shifted away from zero. Bx, Ex and Ey histograms 
are practically indistinguishable. 



3. Comments 



Several comments are in order. 

(i) We performed analogous calculations for CuCu collisions at the top 
RHIC energy and for PbPb collisions at the LHC energy. We found very 
similar qualitative behavior - in fact at t = and x = the following ap- 



proximate scaling holds for (li?^. {\Ex^y\) 



Field Vsr^,,jy ^ 



rrir, 



(3) 



where Ra is the radius of the appropriate nucleus, ^/s is the center of mass 
energy and rUp is a proton mass. In the first approximation, the function 
f{x) is universal. 

(ii) To check the uniformity of the magnetic and electric fields in each 
event, we integrated E and B over the transverse spherical domain Q (de- 
noted in Fig. [T]by a grey circle), i.e.. 



Jn 



(4) 
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and analogously for By^ and Ey. Here r^) is the radius of the domain 
which we varied from r£i = lfmtor/) = 2 fm. As seen from Fig. O for 
6 = 4 fm the ratio R = {\By^Q\) / (|-Ba;,o|) is changing from 1.6 to 2.5 (with 
increasing r^), as compared to 1.14 at x = 0. For 6 = 8 fm i? is changing 
from 3.2 to 5.6, as compared to 1.6 at x = 0. The ratio {\Ey^^\) / {\Ex,n\) is 
always very close to 1, and again the x and y components of E are almost 
identical to the x component of B. 
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Figure 5: The ratio R = (|-By,o|) / (l^aj.ol) as a function of the impact parameter. and 
By are integrated over the transverse domain O of radius ro [fm]. For comparison, we 
also show the ratio at x = 0. 

As seen in Fig. [5], for small b the ratio R weakly depends on the size of 
the domain in contrast to large b, where the td dependence is significant. 
We performed our calculation up to r^) = 2 fm (size of the domain equals 4 
fm), which in peripheral collisions covers almost the whole interaction region. 
In practical applications to the chiral magnetic effect, the fluctuations of the 
magnetic field are important only in the region with a characteristic size 
corresponding to the size of a sphaleron. At least at early stage of a collision, 
where the influence of the magnetic field is significant, this size is expected to 
be of the order of 1 fm (the characteristic size of a sphaleron is less than the 
magnetic screening length, which is (a^T)"^ in the weak coupling limit). This 
suggests, as seen in Fig. 5, that, indeed, (|-Bx|) and {\By\) are of the same 
order and should be treated equally in calculations of the chiral magnetic 
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effect. 

(iii) In an actual experiment, tfie reconstructed reaction plane (the so- 
called event plane) is, in general, tilted with respect to the one defined by the 
geometry of colliding nuclei. Thus a possible experimental measurement that 
is sensitive to the direction of the magnetic (electric) field in each event, will 
reveal the components B'^ {E'^) and By (Ey), that are linear combinations 
of the original ones -B^. (E^) and By {Ey). This mixing provides an addi- 
tional reinforcement to our conclusion that the different field components 
are comparable in magnitude. 

(iv) Since the electric conductivity is significantly larger than the one as- 



sociated with the magnetic field (at least an order of magnitude, see Refs. 13|, 



14| ) it is quite clear that the current created by the electric field Je dominates 
over the one due to the chiral magnetic effect Jb- As seen from Fig. [2] and 
Fig. [3] we expect the ratio to be maximal for central collisions 

and minimal (but still larger than 1) for very peripheral collisions. 

(v) Our results suggest that the strong electric current should not pose 
a significant problem for the correlation observable (cos((/)i + 02 — 2\1/rp)) 



proposed in Ref. [15[ to measure the chiral magnetic effect. Indeed, this 
observable is sensitive to the difference between correlations in-plain and 
out-of-plain and, in the first approximation, Je seems to be the same in x 
and y directions, in contrast to Jb created in peripheral collisions (but not in 
central collision, where we expect jB,x ~ jB,y)- This point, however, should 
be investigated in more elaborated studies, that take into account realistic 
temperature- and density- dependence of the electric conductivity. 

(vi) In principle, by measuring the charge separatior0 in central collisions 
one may perform rough estimate of the electric conductivity a of the mat- 
ter created in heavy ion collisions. Assuming the dominance of the electric 
current in charge separation process in central collisions, and estimating ef- 
fective lifetime teff and spatial extensions of the hot spot, where the electric 
field is almost constant and homogeneous Ees-, we obtain a ~ -, — ^ „ . Here 
Q is the dipole charge and S±_ is the area perpendicular to the dipole axis. 
This problem is currently under our investigation. 

(vii) Our estimates of the magnetic and electric fields are valid only at 
the early stage of the collision. At later stages, the magnetic response from 



the created medium becomes increasingly important [17|, |18| , and may lead 



For instance, in the way proposed in Ref. [If 
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to substantial increase of the magnetic field over the electric one. For quan- 
titative analysis of this effect, more elaborated studies are required 19 . 



4. Summary 

In this Letter we argue that owing to the fluctuating proton positions in 
the relativistic heavy ion collisions, the strength of the magnetic and electric 
fields are comparable (both in the direction parallel and perpendicular to 
the reaction plane) in the event-by-event analysis. Thus, for the observables, 
that are sensitive to the electric and magnetic fields, it is essential to take 
into account event-by-event field fluctuations. 

Our result in combination with the lattice QCD calculations of the elec- 
tric and the chiral magnetic conductivities suggest, that the possible charge 
separation measured separately in-plane and out-of-plane will be dominated 
by the electric current, which in principle allows to measure the electric con- 
ductivity of the hot medium created in heavy ion collisions. However, since 
the electric field is not correlated with reaction plane, the difference between 
correlations in-plain and out-of-plain in peripheral collisions should be insen- 
sitive to the electric field. 
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